RSICC CODE PACKAGE CCC-718

1.         NAME AND TITLE

MCNP-POLIMI EXE :    Monte Carlo N-Particle Transport Code System To Simulate Time-Analysis Quantities.

RELATED DATA LIBRARY

This package includes a test library of cross sections for running the sample problems, but the test library is not suitable for real problems. The user is advised to select available ENDF-based MCNP libraries (60C), or those with the most detailed photon production description for the particular problem. DLC-200 / MCNPDATA (or equivalent) is suitable for use with MCNP-PoliMi. Users who want to receive these data must include MCNPDATA on request, license and Export Control forms.

2.         CONTRIBUTOR

Polytechnic of Milan, Milano, Italy.

3.         CODING LANGUAGE AND COMPUTERS

Fortran 77; Intel-based PCs. (C00718/PC586/01).

4.         NATURE OF PROBLEM SOLVED

MCNP is a general-purpose, continuous-energy, generalized geometry, time-dependent, coupled neutron-photon-electron Monte Carlo transport code system. Based on the Los Alamos National Laboratory code MCNP4C (formerly distributed by RSICC as CCC-700), MCNP-PoliMi was developed to simulate time-analysis quantities. In particular, the code includes the correlation between neutron interaction and the corresponding photon production. Conversely to the technique adopted by standard MCNP, MCNP PoliMi samples secondary photons according to the neutron collision type. A post-processing code, i.e. the Matlab script "postmain", is included and can be tailored to model specific detector characteristics. These features make MCNP-PoliMi a versatile tool to simulate particle interactions and detection processes.

CCC-740/MCNP5 is the most recent MCNP release from LANL and is recommended for general purpose calculations.

5.         METHOD OF SOLUTION

MCNP treats an arbitrary three-dimensional configuration of materials in geometric cells bounded by first- and second-degree surfaces and some special fourth-degree surfaces. For neutrons, all reactions in a particular cross-section evaluation are accounted for. Both free gas and S(alpha, beta) thermal treatments are used. Criticality sources as well as fixed and surface sources are available. For photons, the code takes account of incoherent and coherent scattering with and without electron binding effects, the possibility of fluorescent emission following photoelectric absorption, and absorption in pair production with local emission of annihilation radiation. A very general source and tally structure is available. The tallies have extensive statistical analysis of convergence. Rapid convergence is enabled by a wide variety of variance reduction methods. Energy ranges are 0-60 MeV for neutrons (data generally only available up to 20 MeV) and 1 keV - 1 GeV for photons and electrons.

The MCNP-PoliMi code was developed to simulate each neutron-nucleus interaction as closely as possible. In particular, neutron interaction and photon production are made correlated and correct neutron and photon fission multiplicities were implemented. At each neutron and gamma collision occurring inside user-defined cells, relevant information is recorded, for example reaction type, target nucleus, energy deposited, and position.

The number of photons born in induced fission is determined by sampling from a fission photon multiplicity distribution. The number of photons that can be generated by the algorithm ranges from zero to 23. Photon energy is assigned by independently sampling from an energy distribution taken from CCC-699/MCNP-DSP. In MCNP-PoliMi, a dual-particle source is generated at each induced fission. Neutrons are emitted at the same time with multiplicity and energy distributions given by appropriate distributions from MCNP-DSP. Gamma ray multiplicities and energies have also been implemented on the basis of results found in the literature. This procedure attempts to model the induced fission chains as closely as possible.

6.         RESTRICTIONS OR LIMITATIONS

MCNP-PoliMi requires analog Monte Carlo. The secondary gamma generation performed by MCNP-PoliMi is in general limited by the information present in the ENDF-based MCNP libraries. Do not use delayed neutrons or multigroup cross sections.

7.         TYPICAL RUNNING TIME

Running time varies greatly depending on problem parameters. The cf.inp sample problem describes a Cf-252 point source placed at one meter from a plastic scintillator detector. The detector is cell #1, and the collision data is recorded for this detector. The program will execute this problem in 8.48 minutes on a 2.2 GHz Pentium 4.

8.         COMPUTER HARDWARE REQUIRE­MENTS

MCNP-PoliMi is operable on Windows-based PCs.

9.         COMPUTER SOFTWARE REQUIRE­MENTS

The included Windows executable (file Pcexe\Main.exe) was created with Compaq Visual Fortran Professional Edition 6.6 on a Pentium IV under Windows XP Professional. The developers also tested the code under WINDOWS 98 SE on INTEL Pentium III, under WINDOWS XP Home on AMD Athlon XP, and under Windows2000. The Linux executable was created with GNU GCC version 4.2.2 gfortran on a Opteron cluster under RedHat Enterprise Linux.  The MacOS executables were created using GNU GCC version 4.0.1 gfortran on a MacOS 10.5.8 Operating system.  The MDAS parameter inside the code was set to 16,000,000 (16 Mwords). The plot feature is not included in the distributed executable.

10.        REFERENCES

Adobe Acrobat Reader freeware is available from http://www.adobe.com to read and print the electronic documentation.

a. included documentation in electronic format on the CD in DOC\C718.PDF:

gREADME.txth (June 2011).

E. Padovani and S. A. Pozzi, gMCNP- PoliMi ver. 1.0 Userfs Manual,h Polytechnic of Milan, Italy (November 25, 2002).

b. background information not included in the package :

S. A. Pozzi, E. Padovani, and M. Marseguerra, gMCNP-PoliMi: a Monte Carlo code for correlation measurements,h Nuclear Instruments and Methods in Physics Research (Section A), Volume 513, Issue 3, pages 550-558 (November 2003).

M. Marseguerra, E. Padovani, S. A.Pozzi, M. Da Ros,hPhenomenological simulation of detector response for safeguards experiments,h Nuclear Instruments and Methods in Physics Research, (Section B), Volume 213 pages 289-293 (January 2004).

S. A. Pozzi, E. Padovani, J.K. Mattingly, and J.T. Mihalczo, gMCNP-PoliMi Evaluation of Time Dependent Coincidence between Detectors for Fissile Metal Vs. Oxide Determination,h Institute of Nuclear Materials Management 43rd annual meeting, Orlando, Florida (June 23-27, 2002).

S. A. Pozzi, J.K. Mattingly, J.T. Mihalczo and E. Padovani gValidation of the MCNP-PoliMi code for the simulation of nuclear safeguards experiments on uranium and plutonium metal,h Nuclear Mathematical and Computational Sciences M&C2003, April 6-11, 2003, Gatlinburg, Tennessee, USA.

S. A. Pozzi and J.T. Mihalczo, gMonte Carlo Evaluation of the Improvements in Nuclear Materials Identification System (NMIS) resulting from a DT Neutron Generator,h Institute of Nuclear Materials Management 43rd annual meeting, June 23-27, 2002, Orlando, Florida, USA.

11.        CONTENTS OF CODE PACKAGE

Included are the referenced electronic documents in (10.a), PC executables, and test problems transmitted on CD ROM in 7z format.

12.        DATE OF ABSTRACT

January 2004, June 2011.

KEYWORDS:       COMPLEX GEOMETRY; COUPLED; GAMMA-RAY; MICROCOMPUTER; MONTE CARLO; NEUTRON; SCINTILLATION DETECTOR